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It is known that the ionic conductivity can be obtained by using the diffusion 
constant and the Einstein relation. We derive it here by extracting it from the 
steady electric current which we calculate in three ways, using statistics analysis, an 
entropy method, and an entropy production approach. 
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INTRODUCTION 

Electric conductivity, a, is a coefficient that appears in Ohm's law, 

3 = o-E, 



(1) 



where j is the electric current density and E is the electric field. For ionic conduction, a is 
also called the ionic conductivity. The ionic conductivity can be easily derived by studying 
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FIG. 1: A simple ionic conductor. The white circles are the non-movable lattice, and the black 
circles are the movable interstitial ions that can jump from site to site. For the right graph, there 
is an external electric field pointing upward. The electric field will affect the jump of ions along its 
direction. So one only needs to study how ions jump up and down. 



an ionic conductor that is composed simply of a non-movable crystal lattice and movable 
interstitial ions. The ions can jump from site to site, as shown in Fig. [TJ and external 
electric field E can affect their jumps. If E is chosen pointing upward, one only needs to 
study how an ion jumps up and down. 
Given the ionic jump frequency 



(2) 
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one can get from kinetic theory the diffusion constant [l[ as 

ZW^exp(-^j, (3) 

where v is an effective vibration frequency, e is the height of the barrier which an ion must 
surmount in order to pass to the adjacent site, a is the lattice constant, k B is the Boltzmann 
constant and T is the temperature. The diffusion constant is related to the ion mobility fi 
by the Einstein relation, 

"-iff* (4) 

where q is the interstitial ion charge. Thus one gets the ionic conductivity [l| as 

na 2 q 2 ( e \ 
a = n ^ = l^T UeXP {-k^f)^ (5) 

where n is the ion concentration. One may also obtain the same ionic conductivity by 
calculating the ion drift velocity, 

/ e-Eqa/2\ ( e + Eqa/2\] ( e \ a 2 q ^ 

vd = a r xp 1 — ^f- ) - " exp [ — k^f-)\ w " exp \-ur) w E - (6) 

With the ion drift velocity, one can obtain the electric current density, 

no? q** ( s \ 
j = nqVi = l^T VeXp [-k^f) E ' (?) 

from which one can extract again the same ionic conductivity Eq. 

We shall derive the ionic conductivity again in three different ways, using statistics analy- 
sis, an entropy method, and an entropy production approach. All our derivations are subject 
to three constraints. 

1. N ^> 1: the conductor is large enough to contain a large number (N) of ions. 

2. n ~ 0: the ions are dilute, so each ion jumps independently. 

3. E ~ 0: the external electric field is weak, so the electric current is weak, / ~ 0. 



STATISTICS ANALYSIS 

From Fig. [TJ we understand that when the ion conductor is of small size ( this is the 
same as saying that iV is small ), the electric current will fluctuate. But when the conductor 
is of large size ( i.e., N ^> 1 ), the electric current will become steady. The steady electric 
current is thus actually the most probable electric current, which we will find in this section 
by using statistics analysis. Then we can extract the ionic conductivity. 

For the ionic conductor as shown in Fig. [TJ, we separate the system from the environment. 
The system only consists of all the interstitial ions (and only about how they jump). The 
environment includes everything else. Since we choose n ~ 0, the ions are far from each 
other. Thus, when they jump, they do not affect each other. So we may imagine that all the 
ions jump in a synchronized way and that each jump takes the same time r. This allows us 
to easily study what the jump configuration is and what the corresponding probability is. 
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First, as a start, let the system consist of only one ion. That one ion, in time r, either 
jumps up or jumps down. So we have 
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where 



• f represents the state transition in which the ion jumps up; 

• 4 represents the state transition in which the ion jumps down; 

• Qs is the number of system microscopic states from which the transition can occur; 

N 

• Iv = jV = Yl 1 v ii where V is the conductor volume and v = ±a/r, the effective ion 

i=l 

velocity; 

• Q is the heat exchanged with the environment; 

• ASe is the entropy change of the environment during r; 

• Qe is the number of environment microscopic states compatible with the transition; 

• P is the probability for an ly to occur, and for one ion system P to become the 
Boltzmann distribution; but this P is different from its counterpart in Eq. ([6]) by an 
extra factor 2 which will persist to the end. 

Then, let the system consist of two ions. Eq. ((§]) becomes 



configuration 


tt 


U or it 


n 


Q s 


1 


2 


1 


Iv 


2aq/r 





—2aqjr 


Q 

AS E 
Qe 
P 


2aqE 

laqE 

« ex p T (!f ) 

oc 1 x exp( fcR 9 T ) 





oc 1 
oc 2 x 1 


-2aqE 

2aqE 

« ex P (-« ) 
oclxexp( k ^ T ) 



(9) 



Finally, let the system consist of N ions. Eq. (J5} becomes 
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where k is the number of ions that jump up. The most probable k can be obtained by 
or by 

which becomes 

6\nC k N 2aqE _ n 

Since we choose E ~ 0, we have fc ~ AT/2; and since we choose AT ^> 1, we can apply 
Stirling's formula. So we can write 

InC* «lnC^ - 1 N *> . (14) 
Thus we obtain the most probable fc as 

+ ^ (15) 
^ " 2k B T + 2 • (15j 

Then the most probable Iy is 

Na 2 q 2 E , . 

A; B Tr 

This is also the steady ly. So the steady electric current density is 

. I v na 2 q 2 
J= V= k^fr 

and the ionic conductivity is extracted as 

* = < 18 > 

k B Tr 

The ion jump time relates to the ion jump frequency. Since in this study we only consider 
how an ion jumps up or down, we should have 

(19) 



2z/exp(-^) 



So, Eq. ffl8l) is the same as Eq. (jSJ), except that there is an extra factor 2. 

In our study, we have assumed that all ions are synchronized to jump, and that each 
jump takes the same time interval. In reality, each ion jumps randomly and independently. 
Though the pictures are different, the result is the same. 
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ENTROPY METHOD 

In this section, we shall find the most probable electric current again, but with an entropy 
method. Here the entropy consists of two parts: the environment entropy and the system 
entropy. The environment entropy is obtained from Eq. f TTOT) as 

S E = S Ea + AS E = S m+ {2k -^ a " E , (20) 

where Seo is the initial environment entropy when the ions are not yet ready to jump. Using 
Eq. (fl4j) . the system entropy is 

S s = k B \nQ s = k B \nC k N « S so - BK - 2> , (21) 
where 5*50 is the system equilibrium entropy. The total entropy is 

S = S E + S s . (22) 
Then the most probable k can be obtained by using 

f = ° <»> 

which is equivalent to Eq. ( ITTT) . 

We can also directly calculate the most probable electric current by using 

g = 0. (24) 
To do that, we need write S as a function of I. Using Eq. ( flOi) . 

and using Eq. ff2T]) . the system entropy is rewritten as 

3 > - s » ~ \ W 7 ' (26) 

I — I 

which takes the same form as the entropy of a dilute gas that carries a heat flux [2| or a 
velocity gradient {jjj. Using Eq. (120 \ the environment entropy is rewritten as 

Se = S E0 + (27) 

All the entropies are shown in Fig. |5J Only at the most probable electric current does 
the total entropy have its maximum value. The most probable electric current is also the 
steady electric current. Thus the steady electric current can be obtained by maximizing the 
entropy. It is obtained as 

1 = hdr T K < 28 > 
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FIG. 2: The entropies. The dashed line is the system entropy. The dash-dotted line is the 
environment entropy. The solid line is the total entropy. The total entropy has its maximum value 
at point A. Ia is the most probable electric current, which is also the steady electric current for a 
large enough conductor. 

The corresponding steady electric current density is then 

i = 4 = ^f E , (29) 
J A k B Tr ' v ; 

from which the ionic conductivity Eq. f lXSjl is extracted again. 

In addition, Fig. [2] indicates that the system entropy and the environment entropy 
compete with each other, and their competition determines the value of the steady electric 
current. Entropy always tends to increase. When the environment entropy tends to increase, 
it drives the electric current / to increase too. When the system entropy tends to increase, 
it drives I to decrease in the opposite direction. Only at the steady I do the two entropies 
balance each other, and at the same time, the total entropy has its maximum value. Similar 



discussions can also be found in 17, 18 



ENTROPY PRODUCTION APPROACH 

In this section, we shall see if the steady electric current can be obtained by another 
approach that uses the entropy production. We shall discuss this by analogy with studies 
of the thermal conductivity jl[ and viscosity [3} of a dilute gas. 

For a given electric current /, when it is steady, the entropy production is 

Eh , \ 

°i = -jrl, (30) 

where the subscript I in symbol 07 stands for electric current. When the system relaxes 
from an unsteady state, the entropy production is different. And the shorter the relaxation 
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time, the larger the entropy production. If the relaxation is exponential, and if the minimum 
relaxation time is the jump time r, one may write the electric current that relaxes at the 
maximum rate as 

I fl (t)=Jexp("Y (31) 
Using Eq. (1261) . the system entropy during relaxation is 

1 k B h 2 r 2 2 / 2t 



S s (t) = S so ~ 2jjafr? W I" 7 ' ' (32) 



and the entropy production in the limit t — > is 



<ri 
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t=o ~ Na?q 2±2 ' (33) 



Relations (]HU|) and f[3"3"j) between the entropy production and the current for a steady 
and a relaxing state are shown in Fig. |3j By the steepest entropy ascent (SEA) principle 
we know that an electric current I bigger than I a cannot be steady, because in that 
range the steady electric current does not have the steepest entropy ascent. Thus the steady 
electric current must be in the range [0, I a]- Then the maximum entropy production (MEP) 
principle [jj-16|| indicates that the steady electric current must be I a, because among all 



candidates I a has the maximum entropy production. Thus the steady electric current is 
obtained from I a as 

I= hk^r E ' (34) 
which is the same as Eq. f|28|) . So the same ionic conductivity [18] is reproduced again. 



CONCLUSION 

We derive the ionic conductivity by extracting it from the steady electric current, which 
we obtain by using three different approaches. (I) The first approach is a statistics analysis, 
which is based on calculating the most probable electric current. For a sufficiently large 
ionic conductor, the electric current can be steady, and it is actually the most probable 
electric current, which can be obtained by statistics analysis. (II) The second approach is 
an entropy method. In this approach, we calculate the system entropy and the environment 
entropy, and study their competition. The system entropy is related to the way in which 
each ion moves. The environment entropy is related to the heat generated by the electric 
current. Since both entropies tend to increase, they compete with each other to drive the 
electric current in opposite directions. When they balance each other, the electric current 
becomes steady and the total entropy is maximum. (Ill) The last approach uses entropy 
production principles. In this approach, one needs to assumes that, for a given electric 
current, its relaxation is exponential in time and that the minimum relaxation time is the 
ion jump time. Then, by using two entropy production principles, the SEA principle and 
the MEP principle, one can obtain the steady electric current. 

All three approaches produce the same steady electric current. And so the same ionic 
conductivity is extracted, which is also the same as the kinetic theory result, but with an 
extra factor 2. 
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FIG. 3: Entropy production with respect to electric current. For the solid line, the electric current 
is steady. For the dashed line, the electric current is about to start to relax at the maximum rate. 
The steady electric current is I a- For an electric current I > Ia, it cannot be steady, because in 
this range the entropy production of the steady state is at least lower than that of a relaxing state. 
The steepest entropy ascent (SEA) principle indicates that a state can occur only if it has the 
maximum entropy production. Thus the steady electric current must be in the range < / < I a- 
Subsequently, according to the maximum entropy production (MEP) principle, the steady electric 
current must be I a, because among all candidates it has the maximum entropy production. 
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